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Abstract 

The homotopy analysis method is studied in the present paper. The 
question of convergence of the homotopy analysis method is resolved. It is 
proven that under a special constraint the homotopy analysis method does 
converge to the exact solution of the sought solution of nonlinear ordinary 
or partial differential equations. An optimal value of the convergence con- 
trol parameter is given through the square residual error. An error estimate 
is also provided. Examples, including the Blasius flow, clearly demonstrate 
why and on what interval the corresponding homotopy series generated by 
the homotopy analysis method will converge to the exact solution. 
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1. Introduction 

The search for a better and easy to use tool for the solution of nonlinear equations 
illuminating the nonlinear phenomena of our life keeps continuing. 

A variety of methods therefore were proposed to find approximate solutions. 
One of the most recent popular technique is the homotopy analysis method, which 
is a combination of the classical perturbation technique and homotopy concept 
as used in topology. In the homotopy analysis method, which requires neither a 
small parameter nor a linear term in a differential equation, a homotopy with an 
embedding parameter p £ [0, 1] is constructed. In [TJ a basic idea of homotopy 
analysis method for solving nonlinear differential equations was presented. In this 
method, the solution is considered as the sum of an infinite series, which con- 
verges rapidly to accurate solutions. A numerous nonlinear problems in science, 
finance and engineering were recently treated by the method, see at least [2], [3], 
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0]) 0, |S] j [I], |H],|9] and [TO]. Particularly, a few new solutions of some non- 
linear problems were found by means of the method [TIJ, which were neglected 
by other analytic methods and even by numerical techniques. All of these show 
the potential of the homotopy analysis method for strongly nonlinear problems. 
Recently an optimal homotopy analysis approach for strongly nonlinear differen- 
tial equations were proposed in [12]. Even though a great deal of equations were 
solved using the homotopy analysis method, the question of convergence of the 
method is yet to be answered. 

We in the present paper investigate the homotopy analysis technique from 
a mathematical point of view. The aim is to analyze the method and to show 
that under certain circumstances the homotopy analysis method converges to 
the exact solution desired, without a prior knowledge of the exact solution. An 
optimal value of the convergence control parameter is defined through the square 
residual error concept. Our another emphasis is to address the error estimate of 
the approximate solution. The given theorem is justified exemplifying it by basic 
examples from ordinary and partial differential equations from the literature. 
The presented theory not only gives the convergence, but it also provides the 
information about the interval of convergence of the homotopy series. 

In the rest of the paper, §2] lays the basis of homotopy analysis method. 
A theory is outlined in ^37] for the convergence and error estimate. Illustrative 
examples in §4] are followed by the conclusions in $5} 



2. The Homotopy Analysis Method 

Liao [JJ described the early form of the homotopy analysis method in 1992. The 
essential idea of this method is to introduce a homotopy parameter, say p, which 
varies from to 1 and a nonzero auxiliary parameter so-called the convergence 
control parameter h. At p = 0, the system of equations usually has been re- 
duced to a simplified form which normally admits a rather simple solution. As p 
gradually increases continuously toward 1, the system goes through a sequence of 
deformations, and the solution at each stage is close to that at the previous stage 
of the deformation. Eventually at p — 1, the system takes the original form of 
the equation and the final stage of the deformation gives the desired solution. To 
illustrate the basic ideas of this method, consider the nonlinear boundary value 
problem 

du 

N(u(r)) = 0; r G fl, B(u(r), — ) = 0; r G T (2.1) 

an 

where u(r) defined over the region Q is the function to be solved under the 
boundary constraints in B defined over the boundary V of fl The homotopy 
analysis technique defines a homotopy u(r,p) : R x [0, 1] — > R so that 

H(u,p) = (l-p)[L(u)-L(u )}+phN(u), (2.2) 
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where L is a suitable auxiliary linear operator, u$ is an initial approximation 
of equation (12.11) satisfying exactly the boundary conditions. It is obvious from 
equation ( 12.21) that 

H(u,0) = L(u)-L(u ), H{u,l) = N(u). (2.3) 

As p moves from to 1, u(r,p) moves from Uq(t) to u(r). In topology, this called 
a deformation and L{u) — L{uq) and N{u) are said to be homotopic. Our basic 
assumption is that the solution of equation (12.2ft when equated to zero can be 
expressed as a power series in p 

oo 

u(r,p) = u (r) + pui(r) + p 2 u 2 (r) H = ^u k (r)p k . (2.4) 

The appropriate solutions of the coefficients itfc(r) in (I2.4p can be found from the 
homotopy deformation equations, see [13J. Hence, the approximate solution of 
equation (12.1 ft can be readily obtained as 

oo 

u(r) = \imu(r,p) = uJr). (2.5) 

It was found that the auxiliary parameter h can adjust and control the conver- 
gence region and rate of homotopy series solutions ( 12. 4p . Whenever the series 
( 12. 4 p is known to be convergent, then (I2.5P represents the exact solution of (12. ip . 
as proved in [IB"] . 

3. A convergence Theorem and error estimate 

Using the methodology underlined above, the number of problems treated by the 
homotopy analysis method approaches a couple of hundreds now. However, the 
very basic question of why the series obtained by setting p — 1 in ( 12 .4p should be 
convergent remains unanswered till today. To remedy this issue up to a point, 
we provide the subsequent theorems here. It should be noted that how to find a 
proper convergence control parameter h, or even better, to get a faster convergent 
one, to be used in Theorem 1 will be clarified later in this section. 

Theorem 1. Suppose that A C R be a Banach space donated with a suitable 
norm 1 1| (depending on the physical problem considered), over which the sequence 
Ukit) of ( 12 .4p is defined for a prescribed value of h. Assume also that the initial 
approximation uo(t) remains inside the ball of the solution u(t). Taking r £ R 
be a constant, the following statements hold true: 

(i) If ||f/ c+ i(t)|| < r||t>i;(i)|| for all k, given some < r < 1, then the series 
solution u(t,p) = J2T=o u k(t)p k , defined in (12. 4p converges absolutely at p — 1 to 
( 12. 5 p over the domain of definition of t, 
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(ii) If > T||ffc(t)|| for all k, given some r > 1, then the series solution 

u(t,p) = J2T=o u k(t)p k , defined in (12. 4p diverges at p = 1 over the domain of 
definition of t. 



Proof, (i) In compliance with the ratio test for the power series in p, the 
proof is clear. However, in order to give an estimate to the truncation error of 
homotopy analysis method, we shortly give the whole proof here. If S n (t) denote 
the sequence of partial sum of the series (12.51) . we need to show that S n (t) is a 
Cauchy sequence in A. For this purpose, consider, 

\\S n+1 (t) - S n (t)\\ = \\u n+1 (t)\\ < r\\u n (t)\\ < r 2 |K_i(t)|| < < r n+1 |Mt)||(3.6) 

It should be remarked that owing to (I3.6p . all the approximations produced by 
the homotopy method (I2.2p will lie within the ball of u{t). For every m,n e N, 
n > m, making use of (13. 6p and the triangle inequality successively, we have, 

\\S n {t) - S m (t)\\ = \\(S n (t) - S n -i{t)) + (S n -i{t) - S n - 2 {t)) + ■■■ + {S m+1 (t) - S m (t))\[3.7) 

1 r n~m 

r^\\u Q {t)\\. 

1 — r 

Since < r < 1, we get from (13.71) 

lim \\S n (t) - S m (t)\\ = 0. (3.8) 

n,m— too " v ' v ' " v ' 

Therefore, S n (t) is a Cauchy sequence in the Banach space A, and this implies 
that the series solution (12. 5p is convergent. This completes the proof (i). o 

The proof of (ii) follows from the fact that under the hypothesis supplied in 
(ii), there exist a number I, I > r > 1, so that the interval of convergence of the 
power series (12. 4p is \p\ < l/l < 1, which obviously excludes the case of p = 1. o 



Remark 1. Since the finite number of terms does not affect the convergence, 
Theorem 1 is equally valid if the inequalities stated in (i— ii) are true for sufficiently 
large k's. Thus, it is sufficient to keep track of magnitudes of the ratios rat k 
defined by 



rat k = B^UL , (3 .9) 
and whether they remain less than unity. 

Remark 2. By enforcing the ratio in (i) to hold true in the infinite limit, the 
validity region of t for the series solution can also be constructed, see Example 1 
below. 



Theorem 2. If the series solution defined in (I2.5P is convergent, then it 
converges to an exact solution of the nonlinear problem (12. ip . 
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Proof. The proof can be found in |13j . 



Theorem 3. Assume that the series solution J2^=o u n(t), defined in (12. 5p is 
convergent to the solution u(t) for a prescribed value of h. If the truncated series 
J2n=o u n(t) is used as an approximation to the solution u(t) of problem (12.11) . then 
an upper bound for the error, Euif), is estimated as 

r M+i 

E M {t)<- ||uo(*)ll- (3.10) 

1 — r 



Proof. Making use of the inequality ( 13. Tj) of Theorem 1, we immediately 
obtain 

1 _ r n-M 

\\u(t) - S M (t)\\ < — r M+1 \\u (t)\\, (3.11) 

1 — r 

and taking into account (1 — r n ~ M ) < 1, (13. lip leads to the desired formula (I3.10p . 
This completes the proof, o 

Remark 3. An optimal value of the convergence control parameter h can be 
found by means of the exact square residual error integrated in the whole region 
of interest T, at the order of approximation M, that is, 

. M 

Res(h) = / [NC£u k (r))] 2 dr. (3.12) 

Obviously, the more quickly Res(h) in ( 13 .121) decreases to zero, the faster the 
corresponding homotopy series solution (12. 5 p converges. So, at the given order 
of approximation M, the corresponding optimal value of the convergence control 
parameter h is given by the minimum of Res(h), corresponding to a nonlinear 
algebraic equation of the form 

Therefore, the convergence control parameter obtained via (I3.12H3.1"3"P can be 
supplied into the Theorem 1. However, it is unfortunate that the exact square 
residual error Res(h) defined by (I3.12p needs too much CPU time to calculate 
even if the order of approximation is not very high, and thus is often useless 
in practice. To avoid the time-consuming computation, Liao in [13] suggested to 
investigate the convergence of some special quantities, which often have important 
physical meanings. For example, one can consider the convergence of u'(0) and 
u"(0) of a nonlinear differential equation (12. ip . if they are unknown. It is found 
by the homotopy analysis researchers that there often exists such a region that 
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certain values of h give a convergent series solution of such kind of quantities. 
Besides, such a region can be found, although approximately, by plotting the 
curves of these unknown quantities versus h. These curves are called h-curves or 
curves for convergence- control parameter, which have been successfully applied 
in many nonlinear problems as cited herein. This approach constitutes another 
way of finding proper values of h for the Theorem 1. 

Remark 4. Similar to the constant h— curves idea of Liao, an approximate 
interval of convergence for h can be determined by application of Theorem 1 to 
some certain physical quantities, see Example 3 below. 



4. Illustrative Examples 

To illustrate the validity of the Theorems outlined, we take into account the 
following examples taken from the homotopy analysis studies in the literature. 



Example 1. Consider the first-order nonlinear differential equation [13] 

u' + u 2 = 1, u(0) = 0, (4.14) 

that governs the steady free convection flow over a vertical semi-infinite flat plate 
which is embedded in a fluid saturated porous medium of ambient temperature 
[H] and also the steady-state boundary-layer flows over a permeable stretching 
sheet [15]. In accordance with Theorem 1, choosing uo(t) = 1 — e~ 2t and L = 
4 + 2, the homotopy series solution via the homotopy approach ( 12.21) can be 
straightforwardly constructed. Employing the L 2 norm in R, then Theorem 1 
assures the convergency of (14.141) provided that 



|Mn+1 ^" <1. (4.15) 



\Un(t)\\ 



Table [T] presents evolution of the ratio (14.151) for a variety of the convergence 
control parameter h for the problem (14. 14ft . Additionally, Tabled] gives the error, 
see the last column, defined by 

POO 

err = / \u e (t) - u(t)\dt, (4.16) 



where u e (t) is the solution of (I4.14p calculated numerically. Data displayed in 
Table [T] clearly explains why the homotopy analysis method generates completely 
convergent series solution to the problem ( 14.141) for the chosen parameters h. 

The optimal value of convergence control parameter h is further calculated 
from equations (l3.12H3.13l) only for M = 7 as h = 1.30405, which strongly indi- 
cates the reason of the faster converge of the homotopy series near this value, see 
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h 


M= 1 


M = 10 


M = 20 


M = 30 


M = 40 


M = 50 


err 


2.0 


0.44721 


0.87228 


0.93102 


0.95275 


0.96407 


0.97044 


9.609 x 10 


-5 


1.5 


0.15811 


0.43509 


0.46523 


0.47625 


0.48196 


0.48517 


3.768 x 10- 


-20 


1.3 


0.15811 


0.31877 


0.33354 


0.33881 


0.34152 


0.34304 


8.261 x 10" 


-26 


1.0 


0.31623 


0.45743 


0.47701 


0.48425 


0.48802 


0.49015 


8.384 x 10- 


-18 


0.5 


0.65192 


0.69906 


0.71905 


0.72798 


0.73295 


0.73583 


1.492 x 10' 


-8 



Table 1: The evolution of the ratio (14. 15ft for the equation (14. 14ft . Last column 
corresponds to the error defined by (I4.16p . 



Table [IJ Moreover, the uniform validity region of the homotopy series solution 
(12. 5p can also be analytically evaluated for the special value of h = 1. In this 
particular case, the terms in the homotopy series (12.51) satisfy the ratio 



n+l 



U n (t) 



-2t 



(4.17) 



Hence, regarding (I4.17p . Theorem 1 assures the convergence of the corresponding 
homotopy series (I2.5P for all values of t valid in the interval t > — ^r. 

Example 2. Consider now the second-order nonlinear differential equation 

2u" + u-u 2 = 0, u(0)=0, u(oo) = l, (4.18) 

that governs the steady mixed convection flow past a plane of arbitrary shape 
under the boundary layer and Darcy-Boussinesq approximations [16] . With the 
choices of L = 4^ — 1 and u (t) = 1 — e~*, the homotopy (12. 2 p generates a ho- 
motopy series ( 12. 5 p whose ratios of successive terms are tabulated in Table [5] at 
several order of approximations for some selected values of h. As observed from 
the Table, h = 1.2 generates divergent homotopy solution, and also the conver- 
gence of the homotopy series to the exact solution takes place at a considerably 
slow rate near h = 1 as compared to the smaller values of h. 



h 


M= 1 


M = 10 


M = 20 


M = 30 


M = 40 


M = 50 


err 


1.2 


0.42070 


1.24479 


1.31373 


1.34407 


1.35657 


1.36525 


8828.16 


1.0 


0.24684 


0.88919 


0.94119 


0.96001 


0.96970 


0.97513 


3.444 x 10- 4 


0.8 


0.21484 


0.63937 


0.68278 


0.69872 


0.70700 


0.71167 


1.332 x 10~ 9 


0.6 


0.36425 


0.69932 


0.74524 


0.76240 


0.77137 


0.77643 


1.753 x 10" 7 



Table 2: The ratio of successive terms in the homotopy series corresponding to 
equation (14.181) . Last column corresponds to the error defined by (I4.16p . 
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The optimal value of convergence control parameter h is calculated from equa- 
tions (I3.12ti3.13j) only for M = 7 as h = 0.73258, which explains why the homo- 
topy series (12.51) should converge faster near this value to the solution of (14.181) . 
see Table [2 



Example 3. Consider now the nonlinear partial differential Burger's equation 

u t + uu x = u xx , u(x,0) = 2x, (4.19) 

that has been found to describe various kind of phenomena, such as a mathemat- 
ical model of turbulence and the approximate theory of the flow through a shock 
wave traveling in a viscous fluid [17]. Equation ( I4.19P admits an exact solution 
given by 

2x 

u(x,t) = . (4.20) 

v ' ; 1 + 2* v ; 



To approximate the exact solution (14. 20 p . we choose the auxiliary parameters 

df 



as u (x, t) = 2x and L = Then, the homotopy ( 12. 2 p turns out to be 



Equation (I4.2ip produces the below homotopy series for the solution of (I4.19P 

u(x, t) = 2x - Ahtx + 4ht(-l + h + 2ht)x - Aht(-1 + h + 2htfx (4.22) 
+4ht(-l + h + 2htfx - 4ht(-l + h + 2htfx + ■■■, 

whose convergence to the exact solution (14.201) takes place for the values satisfying 
lim^oo |^j = |1 — h{l + 2t)\ < 1. Thus, in accordance with Theorem 1, this 

holds exactly for the values of — \ < t < together with < h < 2. It is easy 
to demonstrate that the case h = 1 corresponds to the traditional Taylor series 
expansion of the solution (I4.20p . which is only valid in the region — | < t < |. 

Example 4. Consider now the following fourth-order parabolic partial dif- 
ferential equation arising in the study of the transverse vibrations of a uniform 
flexible beam [T8l 



y + z z + x x + y 

Utt + (- 1)U XXXX + (- tyUyyyy + (~ l)u zzzz = (^.23) 

2 cos x 2 cos y 2 cos z 

u(x, y, z, 0) = —u t (x, y, z, 0) = x + y + z — (cosx + cos y + cos z), 

whose exact solution is given by 

u(x,t) = (x + y + z — cosx — cosy — cos z)e~ t . (4.24) 
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To approximate the exact solution (I4.24p . if we choose the auxiliary parameters 

dt 2 ■ 



respectively, u (x, t) = (x + y + z — cosx — cosy — cos z){l — t) and L = J^-, the 



homotopy ( 12. 2ft for h = 1 then generates the subsequent homotopy series 

oo oo / j2n j.2n+l \ 

u(x,t) = ^u n = J2( x + V + z - cosx - cosy - cosz) [ —-^7 - , (ffl ) 

which leads to lim^oo -^y- = 0. So, the homotopy series ( I4.25P converges to 
the exact solution of (I4.23f) for all t. 

Example 5. Let us consider the linear partial differential equation 

u t + u x -2u xxt = 0, u(x,0) = e~ x , (4.26) 

whose exact separable solution is given by 

u(x,t) = e- 8 -*. (4.27) 

To approximate the exact solution (I4.27p . if we choose the auxiliary parameters 
u (x,t) = e~ x and L = respectively, then the homotopy (12.21) turns out to be 

(1 - p)u t (x, t, p) + ph(u t (x, t, p) + u x (x, t, p) - 2u xxt (x, t,p)) = 0, u(x, 0, p) = e"f4.28) 



When h — 1, it appears that the homotopy series solution is convergent only 
for t < 0, that is out of physical interest. On the other hand, for h = —1, the 
successive ratio yields 

, u n+1 (t) . _ \t\ 
u n {t) n + l' 

whose limit gives rise to zero. This explains why the homotopy series solution 
( 12. 5 p in this case represents the real physical solution over t > for h = — 1, as 
also explained in [8]. 

Example 6. As a final example, we consider the classical Blasius flat-plate 
flow problem of fluid mechanics governed by the nonlinear initial- value of third- 
order [19] 

y'" + ^Y = °' y(o) = y'(o) = y '(oo)-i = o, rj e [0,00). (4.29) 

With the transformations 

u . 
V=y t = Ar 7> 

where A is a scaling parameter taken here as 4, system (I4.29P is converted into 
uu" 

u'" + — = 0, u(0) = u'(0) = u'(oo) -1 = 0, t G [0, oo). (4.30) 

9 



Taking L = + with «o(t) = — 1 + 1 + e _t , the values for the convergence 
control parameters h from the homotopy (12. 2 p are shown in Table |3]along with the 
corresponding residual errors. It is seen that the optimal value for the physical 
problem considered is h = —3/2. Figures [fla-d) demonstrate the ratios (13.91) in 
Theorem 1 which were computed from the homotopy series (12.51) for the values 
of h given in Table [3j The convergence is assured further by the values of the 
ratios less than unity, suggesting that limit tends to 0.83, 0.77, 0.76 and 0.78, for 
the values of h shown. 



h 


M = 1 


M = 10 


M = 20 


M = 30 


1.6 


1.258621 


1.201 x 10" 


-3 


2.463 x 10- 


-7 


1.032 x 10- 


-6 


1.5 


1.027450 


1.372 x 10~ 


-4 


3.602 x 10- 


-7 


2.573 x 10- 


-9 


1.4 


0.820771 


1.399 x 10- 


-4 


7.621 x 10- 


-7 


5.454 x 10- 


-9 


1.0 


0.238796 


7.621 x 10" 


-4 


1.511 x 10" 


-5 


4.540 x 10- 


-7 



Table 3: The square residual errors Res(h) for the homotopy series solutions cor- 
responding to equation (I4.30p at several values of convergence control parameters 
h. 



5. Concluding remarks 

In this paper, the homotopy analysis method has been analyzed with an aim to 
investigate the conditions which result in the convergence of the generated homo- 
topy solutions of the nonlinear ordinary and partial differential equations. The 
theorems outlined in the paper have proved that if specific values are assigned to 
the auxiliary parameters in the homotopy analysis method, then the approximate 
homotopy results successfully converge to the exact solution. An optimal value 
approach for the convergence control parameter has also been given. Examples 
have been provided to verify the theory. Via the theorems provided here, not only 
the question of the convergence of the homotopy series is answered, but also the 
region of validity of the space variable ensuring the convergence is determined. 
The traditional flat-plate boundary layer flow problem has been finally treated 
by the convergence theorem. 
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